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Abstract. Small number statistics may heavily affect the 
structure of the broadening function in integrated spec- 
tra of galactic globular cluster centers. As a consequence, 
it is a priori unknown how closely line broadening mea- 
surements gauge the intrinsic velocity dispersions at the 
cores of these stellar systems. We have tackled this gen- 
eral problem by means of Monte Carlo simulations. An 
examination of the mode and the frequency distribution 
of the measured values of the simulations indicates that 
the low value measured for the velocity dispersion of M30 
(Zaggia et al. 1992b) is likely a reliable estimate of the 
velocity dispersion at the center of this cluster. The same 
methodology applied to the case of M15 suggests that the 
steep inward rise of the velocity dispersion found by Pe- 
terson, Seitzer and Cudworth (1989) is real, although less 
pronounced. Large-aperture observations are less sensitive 
to statistical fluctuations, but are unable to detect strong 
variations in the dispersion wich occur within the aperture 
itself. 
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1. Introduction 

Since the pioneering work of Gunn and Griffin (1979), 
much effort has been devoted to investigate the internal 
dynamics of galactic globular clusters (GGCs), both the- 
oretically and observationally. The strategy is to derive 
the fundamental parameters of a GGC, such as mass and 
mass-to-light ratio, by constraining models through the 
radial luminosity and velocity dispersion profiles. The lat- 
ter are measured in two ways: 

i) either from the radial velocities of individual cluster 
stars, or 

ii) from the line broadening of the spectrum integrated 
over some central area. 
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The second approach is more economical than the first 
one, but it suffers from the bias that the integrated spec- 
trum is dominated by few bright stars only, whose light 
overwhelms that of all the others within the spectrograph 
aperture. For this reason it seems important to investigate 
the statistical significance of the measured values a of the 
line broadening and their correspondence to the desired 
estimate of the velocity dispersion {vD^ (see eq. 2). 

The case of M15 is emblematic. Peterson, Seitzer, and 
Cudworth (1989, hereafter PSC) have reported that the 
velocity dispersion at the very center of this core-collapsed 
(King and Djorgovski 1984) GGC is as large as 25 ± 7 
km s~^, a fact implying a steep inward rise of (^2)2. This 
figure results from their interpretation of the measure- 
ments in the framework of the small number statistics; 
in fact, their sampling of (w^)^ in the 3" x 3" central area, 
made using a circular aperture with 2R — 1'.'2, produces 
values of CTcc ranging from 8 to 30 km s~^. (Hereafter the 
subscript cc attached to a indicates a measure of a ob- 
tained with the cross-correlation method; cf. Tonry and 
Davis 1979.) Using a much larger aperture (6" x 6") Du- 
bath, Meylan, and Mayor (1992, hereafter DMM) measure 
directly a lower value CTcc = 14.0 ± 0.3 km s^^. This sec- 
ond result, based on a larger number of stars, is statis- 
tically more significant, but its poorer spatial resolution 
may likely wipe out a rapid increase of (v^)^ in the inner 
1-2 arc seconds, if any. Then the question is: how signifi- 
cant is the peak value of CTcc taken by PSC as the central 
value of (v"^)^ in M15? 

In order to shed some light onto the general problem of 
the statistical significance of the velocity dispersion mea- 
surements from integrated spectra of GGCs, we have in- 
vestigated the effects of the various parameters on sim- 
ulated spectra of the central regions of selected clusters. 
In the next Section we present the methodology, and in 
Sect. 3 we discuss its application to the case of M30, whose 
central velocity dispersion has been measured by Zaggia 
et al. (1992b = Paper I). The case of M15 is discussed in 
Sect. 4. Previous account of this matter has been given by 
Zaggia et al. (1992a). 
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2. Methodology for the simulations 

The width of the broadening function of an integrated 
spectrum represents the luminosity-weighted dispersion 



N 



1/2 



(1) 



of the radial velocities of the stars with luminosity 
ii, falling inside the area viewed by the spectrograph slit. 
(Other effects such as seeing convolution of the spectro- 
graph aperture, convolution of the broadening function by 
stellar and cluster rotation, and presence of close binary 
systems are ignored). Thus a^i may differ from the un- 
weighted radial velocity dispersion for the totality N of 
the stars within the aperture. 
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which is the quantity required in dynamical calculations. 
In fact, the broadening function is dominated by the 
brightest stars. Since they occur in a small number, the 
dispersion of their radial velocities does not necessarily 
represent the average velocity dispersion within the sam- 
pled area. 

We have tested the effect of the small number statis- 
tics and of the star luminosity-weighting in determining 
the value of cr, by means of simulated spectra with a pri- 
ori known {v'^)^] this is the center of the Gaussian dis- 
tribution for {v^)^ produced by poor statistics, i.e. the 
asymptotic value provided by eq. 2 with N oo. Our 
Monte Carlo approach accounts also for the effects of the 
measuring technique used in Paper I, which is based on 
a fit of the peak of the cross-correlation function. Each 
experiment consists of the following steps. 

1. First we generate a random set of stellar luminosity 
values £i with probability distribution given by a lumi- 
nosity function (LF) appropriate to metal poor GGCs 
(Ferraro 1990; see Fig. |). The set is closed when a pre- 
determined total luminosity '^ii = La is reached. 

2. An integrated spectrum is then simulated by sum- 
ming up, for each star of the above set, a same (stan- 
dard) spectrum weighted by the star luminosity £i , and 
shifted in wavelength by an amount corresponding to a 
random radial velocity with Gaussian probability dis- 

tribution of fixed variance (v^ ) oo • 

3. Photon and read-out noise are added to the synthetic 
cumulative spectrum; we decided to tailor these effects 
to the observations of Paper I (global efficiency of the 
3.6 m -I- CASPEC, and characteristics of the CCD used 
there). 

4. Finally, the broadening acc of the synthetic spectrum 
is measured by the same cross-correlation technique 
adopted for the real spectra. 



Input parameters for a set of Monte Carlo simulations are: 

— the total luminosity La = ^£i within the slit aperture 
A; 

— the cut-off value Mum for the LF: stars fainter that 
Miim are ignored in building up the set at step 1. Thus 
Miim is the fixed value which controls the number of 
stars actually concurring to the formation of the sim- 
ulated spectrum; 

— the intrinsic (unweighted) radial velocity dispersion 

{vi)L 

For all the experiments we made use of a high signal-to- 
noise spectrum of HD 122563, the template star of Pa- 
per I, obtained with CASPEC at the ESO 3.6 m telescope 
(resolution of FWHM = 0.156 A, equivalent to 9.0 km s"^ 
at A = 5200 A). We have verified that the results do not 
change by utilizing a lower signal-to-noise spectrum of a 
typical GGC star. We have also tested the use of LFs other 
than the one adopted here (Fig. ||) , finding no appreciable 
difference even with LFs of metal rich GGCs {e.g. NGC 
2808; Ferraro et al. 1990). Finally, simple considerations 
show that no major consequence is expected in the results 
by the fact that stars are assumed to have all the same 
mass and so the same velocity dispersion, provided that 
Miim is faint enough. In fact the difference between the 
mass of a turn-off star and one at 3.5 mag below it, i.e. 
our typical cut-off, is approximately 0.2 M©. With the 
assumption of isothermal velocity distribution, this differ- 
ence corresponds to an increase of ~ 20% of the velocity 
dispersion with respect to that of the turn-off stars. 



Fig. 1. The integrated luminosity function used in the simu- 
lations (solid line) gives the number <E>(M\/) of stars brighter 
than Mv (per magnitude unit). The corresponding luminos- 
ity relative to the total luminosity (in percent) is also shown 
(dashed line). Mum values adopted for M15 and M30 are plot- 
ted as upright dotted lines. 
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Fig. 2. The shaded area in each panel is the frequency distribution of the measured line broadening acc from 1000 Monte Carlo 

experiments with the indicated value of the input velocity dispersion {v^)^ (arrow) and fixed Murn ~ +3.5 mag. The solid line 
is the distribution of the luminosity-weighted d ispGrsion (T ui for the input set of stars (see eq. (1)). The upper left panel reports 
also the histogram of (vr)^ (dotted line) for that particular set of experiments. 



3. The meaning of cTcc- an application to M30 

Here we show how small number statistics and gauging of 
the broadening function operate by reporting on Monte 
Carlo experiments tailored to our observations of M30 
(Paper I). This is a GGC which, by analogy with M15, is 
suspected to possess a central spike in the velocity disper- 
sion (Zaggia et al. 1991). In all simulations La has been 
fixed to —3.5 mag to match the total light (m^ ~ 11.0 
mag) within the centered aperture of 4" 6 x 6" used for 
the spectra of Paper I. Simulations have been made for 
three values of Mum = 1.5, 3.5, and 7.5 mag, and four 

values of the input velocity dispersion {v'^)So — 6, 12, 18, 
and 24 km s~^. 

The results for groups of 1000 simulations are shown in 

the panels of Fig. ||for the four values of {v'^)So (indicated 
by an arrow in the figures), keeping Mum fixed at 3.5 



mag. The shaded area marks the distribution of the output 
values (Tec of the cross-correlation method applied to the 
simulated spectra. We also plot the luminosity-weighted 
velocity dispersion of the input sets of stars (solid line), 
as it results from eq. (1). 

The effect of the small number statistics stands out 
clearly in the upper-left plot of Fig. |2[ where we have 
also plotted the distribution of (w^)^ (dotted line); this 
unweighted equivalent of a^j is calculated for each set of 
stars with eq. (2). The histogram of (v^)^ (truncated in 
its upper part) has a smaller dispersion than those of 
and CTcc. The mean of the distribution of (u^)^ is equal 

to the input value (w^)Sd, and the scatter is due to the 
small number of stars: this implies that the standard de- 
viation of is Poissonian. In fact, with a mean num- 
ber of ^ 74 stars per set (brighter than Mum) and for 
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{v'^)oo — 6 kms^^, the Poissonian standard deviation is 
0.49 km s~^, well in agreement with the standard devia- 
tion of the distribution of {v^)^ derived from our simula- 
tions (0.51 km s^^). 



function. In fact, with a high value of {v^)^ , the cross- 
correlation function become multimodal, with peaks pro- 
duced by bright single stars or groups of stars with radial 
velocity far from the mean. All but one of these peaks 
are ignored by the fitting algorithm, which fits the high- 
est peak only; consequently the measured dispersion is 
systematically lower than the input value. Paradoxically, 
this effect increases with the spectral resolution. 



Fig. 3. Comparison between input ((Ur)^) and output (ctcc) 
values of the velocity dispersion from Monte Carlo simulations 
tailored to our (Paper I) observations of M30. The horizontal 
dotted line represents the measure of CTcc reported in Paper I. 

The luminosity weighting permits some bright stars 
or group of stars to dominate the velocity dispersion. 
So, the distribution of a^, calculated with eq. (1), comes 
out broader than that of {v^)^ ■ For the experiment with 

{Vr)oo = 6 km s~^ , the distribution of a^, has a mean value 

of 5.5 km s^^ and a standard deviation of 1.1 km s~^. This 

behavior is confirmed by the other simulations, and could 

be regarded as the limit to our kinematical measurements: 

in fact, what we measure in a spectrum is a^, which is 

most probably < (f^)^. For instance, in the application 

to the case of M30 we find that the most probable value of 

i_ 

CTu, is lower, by ~ 10%, than the input {v"^)^, and that the 
associated uncertainty (at 1 sigma level) is 20% of ((t«,). 
While the cut-off of the LF moves towards fainter mag- 
nitudes, the intrinsic error lowers. This is the reason why 
we have chosen M^m = 7.5 mag; the integrated light of 
all stars brighter than this limit is 95% of the total light 
encircled by the spectrograph aperture. 

By looking at the successive panels of Fig. ||it is appar- 
ent that the distributions of CTcc and of Gw differ more and 

more as {v"^)^ increases. The distribution of CTcc is also 
increasingly skewed towards the low values. This is due 
to the somewhat unsmooth form of the cross-correlation 



Fig. 4. Mean radial velocity {vr) for each of the PSC spectra 
of M15, scaled to the average radial velocity of 120 single clus- 
ter stars (systemic velocity), plotted against the corresponding 
velocity dispersion (Jcc (large solid dots). Error bars are from 
PSC. The small^dots show the behaviour of 1000 simulations 
made with {v1)'ia ~ 20 km s^^. There is no clear correlation 
between CTcc and {vr). This fact suggests that the largest values 
of (Tec are not likely due to anomalously large radial velocities 
of dominating bright stars. 

The difference between the input and output values 
of the velocity dispersion in these simulations of M30 
can be better appreciated in Fig. ^ where the mode of 
the distribution of dec, together with its half width at 
half maximum, is plotted against the corresponding input 
value of (u^}2 for three values of Mum- It is clear that, 
with the instrumental set up of Paper I, any result giving 
(Tec ^ 10 km s^^ would underestimate, on average , the 
intrinsic velocity dispersion of M30. In Paper I we have 
measured acc = 6.0 ± 0.6 km s~^, well in the range where 
(Jcc — {Vr)^- Note however that, the same figure for acc 
would be (marginally) in accord also with values of (w^}^ 
as high as ~ 25 km s~^ or more. We shall expand the 
discussion on this important point in the next Section. 

It can be of some interest to study the effects of chang- 
ing the size of the integration aperture. We have thus 
investigated the relation between acc and (w^)^ foi' an 
aperture of 1" x 1" (about the same size of that used 
by PSC for M15), which, at the center of M30, encircles 
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a total luminosity La of -0.7 mag (or tua — 13.8). We 
have explored the parameter space using the same values 
of Miim and (w^}^ as in the previous experiments, tak- 
ing care to reject all the experiments in which the total 
light is contributed by one star only. As expected, the de- 
parture of the output from the input velocity dispersion 
increases dramatically. The measured CTcc seems to reach a 
maximum between (w^)^ = 10 and 20 km s~^, then it de- 
creases. Consequently, under these conditions a measure 
of CTcc = 5-^8 km s^^ from a single spectrum would have 
no other meaning but that of a possible lower limit. 



In Fig. ^ we report the results of the simulations for 
the PSC parameter-set. Here {Vr)oo varies from 5 to 40 
km s~^ in steps of 5 km s~^, Mum = 6.8 mag, and La = 
— 1.7 mag {niA = 13.5) as this is the (average) luminosity 
collected by the circular aperture of 1''2 used by PSC to 

sample the velocity dispersion in the core of M15. What 

1^ 

has been just discussed for M30 is still valid: CTcc < {Vr)oo 
for {v'^)So > 10 km s~^ 



Fig. 5. Same as Fig. |, but for M15 as observed by PSC. The 
data refer to the set of simulations with Mum = -1-6.8 mag. The 
horizontal dashed line is the mean of all the PSC values, while 
the dotted line is the mean after rejection of two measurements 
deviating more than 3 sigma. 



Fig. 6. The dotted shaded area in each panel represents the fre- 
quency distribution of acc for 1000 Monte Carlo experiments, 
tailored to the PSC observations of M15, for the indicated value 

of the input velocity dispersion {Vr)So and fixed Mum = -1-6.8 
mag. The solid shaded area is the histogram of the measure- 
ments of M15 velocity dispersion by PSC. 



4. The central cTcc in M15 

The already mentioned claim by PSC of a steep gradient 
of the velocity dispersion in M15, consequence of the high 
value of the central acc > 25 km s~^, seems contradicted 
by the measure of acc = 14.0 km s~'^ reported by DMM. 
We shall note here that the most immediate sign of a 
statistical bias, i.e. the presence of a correlation of the 
velocity dispersion with the corresponding absolute value 
of the excess cluster systemic velocity, is not present in 
the PSC data (Fig. ||). To reconcile the two apparently 
conflicting observations by PSC and DMM, we tailored 
the parameters of our Monte Carlo simulations machinery 
to M15. 



For a better comparison with the PSC observations, in 
the various panels of Fig. ^ we plot the frequency distri- 
butions of acc values resulting from the 1000 simulations 
made for each input value of {v'^)Sa. Superimposed there is 
the histogram of the single measures of acc by PSC (their 
Table 4). It is clear from the panels of Fig. ^ that: 

— Simulations allow to exclude that the central velocity 
dispersion of M15 can be as low as 10 km s~^. 

— There are a chance to measure Ccc > 20 only if 

> 15. 

— The peak of the values measured by PSC is nearly 
coincident with the mode of the distributions of acc 
only when 15 < (v'^)So < 30. 
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The PSC claim of {vD^ > 25 km s^^ at the center of 
M15 rests on the distribution of the results of their eight 
independent measurements of Ccc made around the clus- 
ter center. Six of them lie in the interval 8.4 < acc < 11.8 
km (with a mean value of 10.2 km s~^), but the other 
two are of 23.6 ± 3.1 and 30.0 ± 4.3 km s~^ respectively, 

too large for {v'^)So ^15 km s^^ according to our simu- 
lations. A more quantitative evaluation of the compari- 
son between observations and simulations has been made 
using the Kolmogorov-Smirnov (KS) test. The KS prob- 
ability of the match of PSC data with the Monte Carlo 

n - 

distributions of Fig. ^ is maximum (70%) for {v'^)io = 20 

kms""'^, but it is still large at both {v'^)So = 15 and 25 
km s~^ (63% and 53% respectively). 



meaning. At the same time, we shall see that the aperture 
used by DMM is large enough to wipe out fine kinematical 
features. 

DMM observed M15 with CASPEC through a centered 
aperture of 6" x 6" (as we did for M30). The data plotted 
in Fig. ^ have been computed with the same set of param- 
eters as in Fig. ||but La = —4.5 mag (or ttia = 10.7). We 
see that DMM measure of dec = 14.0 km s^^ is compatible 
with a similar low value for {v'^)^,, but it is consistent with 

values as large as {v'^)Sc = 30 km s~^. This is confirmed 
by the frequency distributions of the CTcc values for fixed 

- rr - 

{Vr)So (Fig. ^. When {v'^)So > 15 km s^\ the measure of 
DMM (arrow) is near the peak of the distribution of ace- 



Fig. 7. Same as Fig. |3|, but for a larger aperture of 6" x 6" used 
for M15 by DMM. Again the data refer to the Mh™ = 4-6.8 
simulation. The horizontal dotted line represents the a^c of 
M15 measured by DMM. 



Fig. 8. Same as Fig. |j, with simulation parameters tailored to 
the observations of M15 by DMM. Their value of acc is marked 
by an arrow. 



In conclusion, our numerical experiments confirm that, 
within R ^ 1"5 from the center of M15, the velocity dis- 
persion is likely 20 km s""'^ (most probable value), in 
fair agreement with the qualitative interpretation given 
by PSC to their observations: ''We conclude that the low 
values for the central dispersion are those dominated by 
one or two stars, and that the higher values are more rep- 
resentative of the dispersion of the background lighf^ (ver- 
batim). Indeed, the multiple sampling of (Jcc is the winning 
strategy to overcome resolution and statistical problems 
at once. Clearly, a single measure of dec through an aper- 
ture as small as the one adopted by PSC would be of no 



5. Conclusions 

From the results of the simulations of the spectra of the 
central regions of GGCs, we can draw the following con- 
clusions: 

1. For M30, our measure of (Tcc — 6.0 ± 0.6 km s~^ (Pa- 
per I) is well in the range where the simulations shows 

that acc ^ {Vr)^- 

2. We confirm quantitatively the results of the discussion 
by PSC about the central velocity dispersion of M15. 
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3. The observations of Ml 5 by DMM are consistent with 

an inward rise of the radial velocity dispersion, al- 
though possibly not as steep as claimed by PSC. 

Finally, we find that a single integrated spectrum through 
a large enough centered aperture is capable of providing 
just a lower value of {v^)^ . The multiple independent sam- 
pling technique pioneered by Peterson et al. (1989) seems 
thus essential to investigate fine kinematical structures of 
the cores of GGCs. 
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